Intrinsically faint quasars: evidence for meV axion dark matter in the Universe 
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Growing amount of observations indicate presence of intrinsically faint quasar subgroup (a few % 
of known quasars) with noncosmological quantized redshift. Here we find an analytical solution of 
Einstein equations describing bubbles made from axions with periodic interaction potential. Such 
particles are currently considered as one of the leading dark matter candidate. The bubble interior 
possesses equal gravitational redshift which can have any value between zero and infinity. Quantum 
pressure supports the bubble against collapse and yields states stable on the scale more then hun- 
dreds million years. Our results explain the observed quantization of quasar redshift and suggest 
that intrinsically faint point-like quasars associated with nearby galaxies are axionic bubbles with 
masses 10 8 -10 9 M© and radii 10 3 -10 4 7? Q . They are born in active galaxies and ejected into sur- 
rounding space. Properties of such quasars unambiguously indicate presence of axion dark matter 
in the Universe and yield the axion mass m m 1 meV, which fits in the open axion mass window 
constrained by astrophysical and cosmological arguments. 



Based on observations, Karlsson pj has noted division 
of quasars (QSOs) into two groups with different red- 
shift properties and concluded the following. If we select 
QSOs associated with most nearby (distance d < 50 — 100 
Mpc), galaxies then their redshift is close to certain val- 
ues (quantized), as shown in Fig. 3 below. Meanwhile, 
in QSO samples associated with distant galaxies no pe- 
riodicity in intrinsic redshift is observed. Such a division 
is supported by later studies of QSOs associated with 
most nearby galaxies where the quantization was con- 
firmed 0,13 and distant (0.01 < z ga i < 0.3) galaxies for 
which absence of any periodicity was claimed 0. The 
observations suggest existence of intrinsically faint (op- 
tical luminosity L = 10 5 — 10 7 L Q ) QSO sub group with 
quantized noncosmological redshift. Being intrinsically 
faint, such objects are not detected from large distances 
(which yields disappearance of redshift quantization in 
distant QSO samples) and constitute only a few % of the 
known QSO population. Observations indicate that such 
quasars are ejected from nearby active galaxies or in the 
process of ejection from the galactic nucleus 0,0- 

Here we show that bubbles of dark matter with pe- 
riodic interaction potential, masses about 10 8 — 10 9 M Q 
and radii 10 3 — lO 4 i?0 can explain the intrinsically faint 
quasars. The bubble is supported against collapse by 
quantum pressure and decays on a time scale more than 
hundreds million years. Hypothetical axions, one of the 
leading dark matter candidate, fit well into this picture 
and can account for the redshift quantization. Usual 
baryonic matter falls into the bubble interior, heated 
by the release of the gravitational energy and produce 
electromagnetic radiation that freely propagate into sur- 
rounding space. 

In this Letter we study massive real scalar field tp with 
periodic interaction potential 



V(tp) = V [l-co S (tp/f)], 



(1) 



where Vq > 0. This potential is quite general and de- 
rived in quantum filed theory in connection with pseudo 
Nambu-Goldstone bosons (PNGBs) 0- In all such mod- 



els, the key ingredients are the scales of global symme- 
try breaking / and explicit symmetry breaking (Vo) 1 ^ 4 - 
One of the examples of a light hypothetical PNGB is 
the axion which possess extraordinarily feeble couplings 
to matter and radiation and is well-motivated dark mat- 
ter candidate 0. If the axion exists, astrophysical and 
cosmological arguments constrain its mass to be in the 
range of m — 10~ 6 — 3 x 10~ 3 eV and the global 
symmetry-breaking scale to lie in a window / « 10 7 GeV 
x0.62/m(eV) = 2 x 10 9 - 6 x 10 12 GeV @. 
We consider spherically symmetric system with metric 



ds 2 = ~N 2 dt 2 + g 2 dr 2 



*dn 2 



(2) 



where g, the radial metric, and A, the lapse, are functions 
of t and r with r being the circumferential radius. We 
introduce dimcnsionlcss coordinates and define the unit 
of distance, time and tp as 
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where c is the speed of light, G is the gravitational con- 
stant, m = y/Vo/f is the particle mass. In dimensionless 
units the static Klein-Gordon and Einstein equations de- 
scribing the self-gravitating field tp and the metric are 
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with boundary conditions 

g(0) = g(oo) =A(oo) = l, 
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ff '(0) = iV'(O) = ^'(0) = 0, V(<p(oo)) = 0, 



where prime denotes d/dr, 



V 



or 



cos(ay)], a 



'pi 



(J) 



is the dimensionless potential and the coupling parame- 
ter respectively, m p \ = y/Tic/G = 1.2 x 10 19 GeV is the 
Planck mass. The interaction potential has degenerate 
minima at ip — 2im/a, where n is an integer number. 
Here we show that in the limit of strong nonlinearity, 
a ^> 1, Eqs. |@J-© have an approximate static solu- 
tion that describes a spherical bubble with surface width 
much smaller then its radius R. The bubble surface is 
an interface between two degenerate vacuum states with 
ip = 2wn/a (r < R) and <p — (r > R). 

Outside the bubble Eqs. (0}-© lead to the known 
Schwazschild solution 
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iV 2 = l-^, (8) 
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where M is the bubble mass in units of m^j/m. 

Let us assume that JJ»a» 1. Then, near the surface 
one can omit terms with X/r in Eqs. and take 

r ~ R, we obtain 



w" dV 
2R<p'V + - — = 0, 
9 dip 



N'=™ (^' 2 - 2g*V) , 



g ' = i^^ 2 + 2g 2 V) 



(9) 



(10) 



(11) 



Eqs. 101- dJ can be solved analytically. Their first inte- 
gral is 



N = const, ip' 2 = 2g 2 V, g' = Rgtp K 



(12) 



We assume tp(0) — 2im/ot, where n = 1,2,3, .. . is the 
number of kinks at the bubble surface, and ip(r) mono- 
tonically decreases with r. Eqs. (|12fl yield 
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For V{tp) given by Eq. J7|) the final solution is 
4i? 
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tp G [27r(n - m + l)/a,27r(n - m)/a], (14) 
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FIG. 1: Scalar field ip as a function of distance r to the bubble 
center for bubbles with equal radius and different quantum 
numbers n = 1, 2, 3. The unit of length is %/mc. Note, 
we plot the field <p only in the vicinity of the bubble surface 
where it undergoes variation. 



where R m is a position of the 771th kink, m — 1,2, ... ,n. 
When the coordinate r passes through the point R m 
the scalar field ip(r) changes from 27r(n — m + l)/a to 
27r(n — m)/a (see Fig. 1). Eq. i|13|) yields the following 
expression for g as a function of (p inside the bubble: 



1 4i? r 
- = 1 - -H2m- 
9 " 



1 + cos(a^/2)]. 



(15) 



Outside the bubble ip — 0, m = 77 and X/g = 1 — 
8nR/a 2 . The solution is valid if X/g > 0, that is 
R < i?max = a 2 /in. Match of the inner solution il 1 -il) 
with the Schwazschild solution © determines the mass- 
radius relation 



M = AirnuR 2 - 8ir 2 n 2 u 2 R 3 , 



(16) 



where u is the surface energy density given by an integral 
over one potential period u = J \J 2V dip / An . For the 
cosine potential u = 2/ na 2 . 

Redshift of the bubble interior z = X/N — 1 can be 
found by matching the inner N — const and the outer 
(jHJ) solutions: 
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The internal redshift monotonically increases from zero 
to infinity when the bubble radius R changes from zero to 
Rmax- Fig. 2 shows the redshift of space as a function of 
the distance r to the bubble center. The redshift is con- 
stant in the bubble interior and monotonically decreases 
outside the bubble. 

Let us make rescaling M — > M/Attu, R — > R/Attu, then 
Eqs. (TBJ}, 113 yield 



M = nR z - n 2 R 6 /2, 



X 



1 -nR 



X. 



(18) 
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FIG. 2: Redshift z of space as a function of distance r to the 
bubble center for bubbles shown in Fig. 1. 



For a given M the redshift depends on the integer number 
n, which implies the redshift is quantized. 

In early samples of QSOs associated with nearby spiral 
galaxies, Karlsson showed that the redshift distribution 
has a periodicity log(l + z n+1 ) — log(l + z n ) — 0.089, 
where n = 0, 1, 2, . . . and z = 0.061 0. It has been 
later confirmed by other groups 0, El • In a recent pa- 
per, Burbidge and Napier |3j tested for the occurrence of 
this periodicity in new QSO samples and found it to be 
present at a high confidence level. The peaks were found 
at z w 0.30, 0.60, 0.96, 1.41 and 1.96 in agreement with 
Karlsson's empirical formula. The formula also includes 
the peak at zo — 0.061, however, this peak does not occur 
for quasars, but for morphologically related objects. 

The redshift periodicity is observed only in QSO sam- 
ples satisfying certain selection criteria, in particular, the 
galaxies which are assumed to be paired to the QSOs 
must be most nearby spirals Q, Il2| . This implies that 
redshift quantization is a property of intrinsically faint 
QSOs which are not detected from large distances. 

It is naturally to assume that QSOs born in the same 
type of galaxies have approximately equal masses because 
their formation mechanism must be similar. Such phe- 
nomenon is well known for type la supernovae or neu- 
tron stars: practically all measured neutron star masses 
cluster around the value of 1.4M Q with only a few per- 
cent deviation [l3T | . If dark matter bubbles are born with 
equal masses then, according to Eqs. H18|) . their red- 
shift must be quantized. For M = 0.0601 (in dimension 
units M = 0.00752a 2 m 2 )1 /m) Eqs. {TBJ nav <2 solutions for 
n = 1, 2, . . . , 8, they are given in Table 1. 

In Fig. 3 we plot the most recent histogram of the 
redshift distribution from Ref. |T3| in which five peaks 
are clearly seen. The solid lines show the redshifts from 
our Table 1, they match well the observed peaks. The 
agreement is remarkable because the theory has only one 
free parameter, the bubble mass M. Such coincidence 
strongly suggests that the point-like quasars associated 
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R, in a. 2 h/mc 
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0.0329 


0.357 
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0.0159 
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0.0153 
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8 


0.0151 


26.6 



TABLE I: Redshift of the bubble interior z and its radius R 
for M = 0.00752a 2 m^/m and different kink numbers n. 
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FIG. 3: Histogram of the redshift distribution of QSOs close 
to bright nearby active spiral galaxies or multiple QSOs with 
small angular separation from Ref. 12]. The solid lines rep- 
resent position of the peaks from Table 1. 



with nearby galaxies are dark matter bubbles composed 
of scalar particles with periodic interaction potential. 
One should mention an alternative possibility of quasar 
evolution. Bubbles can be originally born with the same 
mass and number of kinks n = 5 that corresponds to 
the 5th peak. During evolution the kinks tunnel to the 
bubble center and quasars sequentially decay into states 
with smaller n but the same mass. 

For axions with m = 0.1 — 3 meV and / = 2 x 10 9 — 
6 x 10 10 GeV Eq. © yields a = 5.6 x 10 7 - 1.7 x 10 9 . 
Hence, an axion bubble with the internal redshift z — 
0.36 and n = 1 would have the mass M — 3 x 10 7 — 
10 9 M Q and the radius R = 3 x 10 2 - 10 4 i? Q . Such radius 
range agrees with the size of the emission region expected 
for the intrinsically faint QSOs. Indeed, for Seyfert 1 
galaxies the size of the broad-line region is R ~ 10 — 100 
light days. The luminosity L of the QSOs is 5-6 orders 
smaller. Based on the empirical relation for Seyfert 1 
galaxies R oc L 1 / 2 jlj], we obtain for the quasars R ~ 
10 4 i?„. 
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Now we discuss the bubble life time. Under the in- 
fluence of surface tension and gravitational attraction an 
initially static bubble starts to collapse. In the thin-wall 
approximation the initial acceleration is given by |l5j 



R = 



2N 3 



N 2 M 



(l + iV)i? 2 ' 



(19) 



where N = yl — 2M/R. For one kink thin-wall con- 
tracting bubble the conserved mass is |lq 



M 



AiruR 2 



V 7 ! - {dR/di 



8tt 2 u 2 R 3 



(20) 



where r is the interior coordinate time. Based on Eqs. 
(|19H , (|20[) one can expect a continuous contraction of the 
bubble to the origin on an astronomically short time scale 
R/c <C lyr. However, so far we treated the scalar field 
as classical. Quantum corrections suppress the collapse 
and result in appearance of long-lived bubbles, stable on 
a scale more then hundreds million years. To include 
quantum effects it has been suggested that the expres- 
sion H20|l be interpreted as the canonical hamiltonian of 
the bubble at the quantum level E3, O- The bub- 
ble wave function *B(R) satisfies the following stationary 
quantum mechanical equation in one dimension (U = 1) 

m 



(E + 8ir 2 u 2 R 3 ) 



d 2 

dR 2 



167rVi? 4 



V(R) = 0. (21) 



This equation possesses stationary solutions that are not 
possible in the classical model. Bubbles of non-negligible 
redshift correspond to highly excited stationary states 
for which the energy spectrum can be treated as quasi- 
continuous. At the quantum level the collapse is pre- 
vented by quantum pressure that balances the surface 
tension and gravitational attraction producing station- 
ary configurations. 

Let us estimate the decay time of an excited station- 
ary state of the quantum bubble. The decay occurs by 
means of scalar particle emission. We estimate the de- 
cay time using the Bohr correspondence principle as the 
time of energy loss by the classical bubble with the ra- 
dius R(t) oscillating between the turning points R(t) = R 
and R(t) — 0, where R is determined by Eq. (|16|) . In 
the quantum picture, however, there are no such oscil- 
lations. The probability of creation a particle with the 



energy mc 2 by a moving bubble surface is governed by the 
Boltzmann factor exp(— mc 2 /T c g ), where T c ff = d/c is 
the effective temperature and a is the acceleration of the 
surface [l^. For the bubble a « c 2 /R(t) and the Boltz- 
mann factor reduces to exp(— R(t)/l) where I = h/mc is 
the surface width. Hence, emission of scalar particles is 
exponentially suppressed apart from small regions where 
R(t) < I. As a result, during one period of oscillation, 
t c ~ R/c, the energy loss is AE ~ (l/R)E, which yields 
the bubble life time 



R R 2 , . 

*~7 <c "' (22) 

For an axionic bubble with R > 10 2 -R© and I < 0.7 cm 

Eq. IL'l'l) yields t > 10 8 yrs which is the time we need to 

account for the phenomenon of quasars. 

Properties of the intrinsically faint point-like QSOs, 
combined with equations for the bubble mass M — 
0.00752a 2 m 2 1 /m = 2.94m(eV) x lO n M and the radius 
R = 0m29a 2 h/mc = 2.73m(eV) x 10 6 i? o , allow us to 
determine the axion mass m. The quasar luminosity sug- 
gests that the bubble radius is larger then lO 3 i?0 which 
yields m > 0.4 meV and M > 1O 8 M . From the other 
hand, the quasar ejection from active galaxies implies 
that the bubble mass M must be much smaller then the 
galactic mass. It is reasonable to constrain M < 10 9 M Q 
which leads to m < 3 meV and R < 10 4 i?Q. We con- 
clude, the axion mass is m = 0.4 — 3 meV. This value fits 
in the open window for the axion mass constrained by as- 
trophysical and cosmological arguments Q , which unam- 
biguously points towards the axionic nature of dark mat- 
ter composing the intrinsically faint point-like quasars. 
Current cavity search experiments in Livermore |20| and 
Kyoto University @ are looking for the axion in the mass 
range 1 — 10 fj,eV which deviates by two orders of mag- 
nitude from our result. Probably now, when the axion 
mass is established from quasar observations, the axion 
has a better chance to be discovered. 

We mention that data on central "black hole" masses 
in small companion galaxies allow us to determine the 
axion mass more accurately and yield m = 1.0 — 1.9 
meV. We will discuss this elsewhere. Moreover, observa- 
tions show that apart from the intrinsically faint point- 
like objects considered here there is a subgroup of bright 
quasars which probably also possess noncosmological red- 
shift. Tachyons, another dark matter candidate, can ex- 
plain their nature. We discuss this in a detail paper [2l| . 
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